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LOCAL INVERSE ESTIMATES FOR 
NON-LOCAL BOUNDARY INTEGRAL OPERATORS 

M. AURADA, M. FEISCHL, T. FUHRER, M. KARKULIK, J. M. MELENK, AND D. PRAETORIUS 


Abstract. We prove local inverse-type estimates for the four non-local boundary integral 
operators associated with the Laplace operator on a bounded Lipschitz domain Q in for 
d > 2 with piecewise smooth boundary. For piecewise polynomial ansatz spaces and d S 
{2, 3}, the inverse estimates are explicit in both the local mesh width and the approximation 
order. An application to efficiency estimates in a posteriori error estimation in boundary 
element methods is given. 


1. Introduction 

Inverse estimates are general tools for the numerical analysis of discretizations of partial 
differential equations (PDEs). They provide a means to bound a stronger (semi-) norm of 
a discrete function by a weaker norm up to some negative power of the mesh width. For 
example, in the context of finite element methods, it is textbook knowledge that 

||^V 14 ||l 2 (q) < C ||14||l2(o) for all continuous T^-piecewise polynomials I 4 . (1.1) 

The constant C > 0 depends only on the shape regularity of the underlying triangulation 
7h of C and the polynomial degree of I 4 . Here, h G L°^{Q) is the local mesh width 
function defined by h\T '■= diam(T) for T G Th- Inverse estimates have also been derived 
for fractional-order Sobolev spaces [GHS05, DFG’''04]. The usual proof of inverse estimates 
like (1.1) relies on scaling arguments, i.e., the powers of h arise by elementwise, i.e., local 
considerations and transformations to reference configurations. 

In the present work we consider the four classical boundary integral operators (BIOs) 
associated with the Laplacian, e.g., the 3D simple-layer integral operator 

If 1 

= — 1 -r My) dy for x G dVL. (1.2) 

47r Jan \x - y\ 

Here, D C M'^, d > 2, is a bounded Lipschitz domain with piecewise G^-boundary dVL. Let 
r C dVL be a relatively open subset of the boundary dVL. Our main result for QJ and d G {2, 3} 
reads, simplified, 

||hi/2(p + i)-iVr5J$,|U2(r) <G||<l>;,||5_V2(r) (1.3) 

for all 7h-piecewise polynomials of degree p G No, where Vr(-) denotes the surface 
gradient. The constant G > 0 depends only on the shape regularity of the underlying trian¬ 
gulation Th of r. In typical settings, QJ is an isomorphism between H~^^‘^{T) and 
so that we observe that (1.3) is in fact an inverse estimate for the finite dimensional space 
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{Q3$/i : a Tft-piecewise polynomial of degree p G Mq} for the weighted iJ^-seniinorm and 

the natnral hT^/^-norm. Inverse estimates of the form (1.3) will be shown for all fonr BIOs 
associated with the Laplacian and discrete spaces with spatially varying polynomial degree, 
cf. Corollary 3.2. In fact, in Theorem 3.1 we will show more general resnlts of the type 

||0||p-i/2(r) + ||wh0||L2(r) for all 0 e L2(r), (1.4) 

L°°(r) 

where Wh is a fairly general weight fnnction. The correct choice of the weight fnnction Wh 
and an inverse estimate from [GHS05, KMR14] for the weighted L^-norm allows one to infer 
(1.3) from (1.4). 

Applications. The inverse-type estimate (1.3) arises natnrally in adaptive BEM (bonndary el¬ 
ement method) when one tries to transfer the convergence and qnasi-optimality analysis from 
adaptive FEM [CKNS08, Ste07] to adaptive BEM [FKMP13, Ganl3]. Indeed, the present 
resnlts allow ns to prove qnasi-optimality of adaptive BEM for piecewise smooth geometries 
and higher (hxed) order discretizations; we refer to [FEK+M] and [FFK’''15], where this is 
worked ont in detail for weakly singnlar and hypersingnlar integral eqnations, respectively. 
While the inverse estimate (1.3) featnres prominently in the analysis of qnasi-optimality of 
adaptive BEM for symmetric problems, it is also a key ingredient for plain convergence in 
non-symmetric problems snch as FEM-BEM conplings. We refer to [AFF+13a] and the ear¬ 
lier preprint [AFF’''12] of the present work for a convergence proof of the adaptive conpling 
of FEM and BEM. 

A fnrther application of estimate (1.4) concerns the efficiency of weighted residnal error 
estimators for both weakly singnlar and hypersingnlar integral eqnations [Car97, CMSOl, 
CMPS04]. To £x ideas, consider the weakly singnlar case and snppose that 0 G T^(P) solves 
QJ0 = / for some given / G 77^ (P). Let be the Galerkin approximation of 0, where the 
ansatz space consists of 7/i-piecewise polynomial of fixed degree p G Nq. While reliability 

Cy 11^ - ■S>/.llg-.«(r) < iKm -■= l|ft‘''Wr(/ - OtJIU^ir) (1.8) 

is well-known (at least for polyhedral domains 72), the converse efficiency estimate remained 
open. As a conseqnence of (1.4), we will see in Corollary 3.4 that 

C~gr]h,m < ||h^'^^(0 - <I)ft)||L2(r), (1.6) 

which expresses efficiency of the weighted residnal error estimator with respect to the slightly 
stronger norm ||h^/^(0 — <l>/i)||i 2 (r) > ||0 — We refer to Corollary 3.7 for the case 

of the hypersingnlar operator. 

These efficiency bonnds are specific instances of new stability estimates for the BIOs in 
locally weighted L^-norms detailed in Corollaries 3.3 and 3.6. 

Novelty. The discrete ineqnality (1.3) was first shown independently in [FKMP13] and 
[Ganl3], however, nnder some restrictions. The work [FKMP13] considers only lowest- 
order polynomials, i.e., 77-piecewise constants, bnt works for polyhedral bonndaries P. The 
work [Ganl3] proves (1.3) for arbitrary 7h-piecewise polynomials, bnt its wavelet-based anal¬ 
ysis is restricted to C^’^-bonndaries P and the constant C > 0 depends on the polynomial 
degree. Onr proof of (1.4) generalizes the works [FKMP13, Ganl3] in the following ways: 
1) we generalize the analysis of [FKMP13] for the simple-layer operator TT to all fonr BIOs 
associated with the Laplacian (i.e., the donble-layer operator its adjoint R', and the 

2 


||M;ftVrT10||L2(r) 


< 


Wh 


£ 1/2 







hypersingular operator QU); 2) we extend our previous analysis from polyhedral domains 
to piecewise smooth geometries; 3) we lift the restriction to hxed-order polynomial ansatz 
space and permit very general ansatz spaces; 4) for ansatz spaces of piecewise polynomials of 
arbitrary order, we make the dependence on the polynomial degree in the inverse estimates 
explicit. 

The technical difficulty in the proof of (1.4) and (1.3) lies in the non-locality of the 
boundary integral operator TT, which precludes simple elementwise considerations. We cope 
with the non-locality of the BlOs by splitting them into near-held and far-held contributions, 
each requiring diherent tools. The analysis of the near-held part relies on local arguments 
and stability properties of the BlOs. For the far-held part, the key observation is that the 
BlOs are derived from two volume potentials, namely, the simple-layer potential QJ and the 
double-layer potential ^ by taking appropriate traces. Since these potentials solve elliptic 
equations, “interior regularity” estimates are available for them and trace inequalities imply 
corresponding estimates for the BlOs. Section 4 proves the relevant estimates for the simple- 
layer potential 03, whereas Section 5 is concerned with the double-layer potential .ft. The 
hnal Section 6 then combines these results to give the proof of Theorem 3.1. 

Although the present paper considers only the four BlOs associated with the Laplacian, 
the scope is wider. As just mentioned, the key tool are interior estimates for potentials; such 
estimates are available for many elliptic equations, for example, the Lame system, so that 
we expect that corresponding results can be proved as well for BlOs associated with these 
problems. 

General notation. We close the introduction by stating that | ■ | denotes, depending on the 
context, the absolute value of a real number, the Euclidean norm of a vector in M'^, the 
Lebesgue measure of a subset of or or the (d — l)-dimensional surface measure of a 
subset of do.. The notation a < b abbreviates a < C ■ b for some constant O > 0 which will 
be clear from the context, and we write a ~ 6 to abbreviate a ^b ^ a. 

We write Br{x) := { 2 ; G : \x — z\ < r} for the closed ball with radius r and center x. 

2. Spaces, Operators, and Meshes 

2.1. Sobolev spaces, is a bounded Lipschitz domain in d > 2, with piecewise C^- 
boundary dQ and corresponding exterior domain := \ 12. The exterior unit normal 

vector held on dQ is denoted by u. Throughout, we will assume that T C dfl is a non-empty, 
relatively open set that stems from a Lipschitz dissection dQ = ru9ru(cA2\r) as described 
in [McLOO, pp. 99]. Note that T = dfl is valid. 

The non-negative order Sobolev spaces for s G { — 1/2, 0,1/2} are dehned as 

in [McLOO, pp. 99] by use of Bessel potentials on and lifting via the bi-Lipschitz maps 

that describe dil. We also need the spaces and In accordance with 

[McLOO], these are dehned as follows: 

:= {w|r: v G (2.1) 

:= { v : Eo,rV G (2.2) 

where Eop denotes the operator that extends a function dehned on T to a function on 
dQ by zero. These spaces are endowed with their natural norms, i.e., the quotient norm 
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||w||j:^i/ 2 +s(r) := inf{||y||^i/ 2 +.(gf^): V\r = v} and ||n||jji/ 2 +.(r) := ||^o,rw||i/i/ 2 +=(aQ)- Owing to 
the assumption that dQ = F U SF U {dQ \ F) is a Lipschitz dissection, we have the following 
facts, stated here without proof: 

Facts 2.1. (i) For s = 1/2, we have the norm equivalences ~ ||w||| 2 (gQ) + 

II VrM||i2(gQ) and ||m||~ + II VrM|li2(r), where Vr is the (weak) surface 

gradient. 

(a) For s = 0, the norms ||M||//i/ 2 (aQ) and ||m||^i/ 2 (p) can equivalently he described by the 
Aronstein-Slobodeckii norms of u and -Eopw (cf. [McLOO, (3.18)] for the definition of 
the Aronstein-Slobodeckii norm). 

(Hi) For s = 0, the spaces H^/‘^{dVt) and are obtained from interpolating between 

the cases s = —1/2 /he., L‘^{dQ) or F^(F )) and s = 1/2 /he., or H^{T)) using 

the K-method (cf., e.g., [McLOO, Thni. B.ll] for the case of H^{dVt)). 

Negative order Sobolev spaces are dehned by duality, namely, for s G { — 1/2, 0,1/2}, 

H-^/\dn) := 

^-(i/ 2+^)(F) := i7i/2+^(F)', and (2.3) 

if-{i/ 2 +^)(F) := 

where duality pairings (•, •) are understood to extend the standard L^-scalar product on dVL 
or F. We observe the continuous inclusions 

^±(1/2+^) (F) C h 7 ±(V 2 +h(r) as well as 

We also note that for G L2(F) the zero extension i?o,r'0 satisfies i?o,r'0 £ 
with 

llV'llH-i/ 2 (r) = ll-E'o,rV'llH-i/ 2 (an)- (2-4) 

We denote by 7 ™*(-) the interior trace operator, i.e., is the restriction of a function 
u G H^{Vt) to the boundary dVt. With H\{Vt) := G H^{Vt) : -Aw G ^2(11) j, the interior 
conormal derivative operator 7 ™* ; Fl^iVt) —)■ Fl~^^‘^{dVL) is dehned by the hrst Green’s 
formula, viz., 

( 7 ”*’m , n)aQ = (Vu , Vv)n - {-Au , v)n for all v G (2.5) 

Remark 2.2. The operator 7 )“*(-) generalizes the classical normal derivative operator: if 
u G FI\{VL) is sufficiently smooth near a boundary point xq, then can he represented 

near xq by a function given by the pointwise defined normal derivative d^u. ■ 

The exterior trace 7 “* and the exterior conormal derivative operator 7 ®”^* are dehned 
analogously to their interior counterparts. To that end, we hx a bounded Lipschitz domain 
U with fl G U. The exterior trace operator 7 “* : F[^{U \ G) —)■ F[^/‘^{dVt) is dehned by 
restricting to dVL, and the exterior conormal derivative is characterized by , n)ao = 

(Vm , yv)u^ — {—Au , for all n G iF^(17 \ G) with = 0 on dU. 

For a function u that admits both conormal derivatives or both traces, we dehne the jumps 
[ 7 im] := 7 “W — 7 ™W and [u] = — 7 ™*m, respectively. 
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2.2. Boundary integral operators. We briefly introduce the pertinent boundary integral 
operators and refer to the monographs [McLOO, HW08, SSll] for further details and proofs. 
Green’s function for the Laplace operator is given by 


G{x,y) 


“ 1 ^ log lx- 2 / 1 , ford = 2, 
+ |g^ |x - ford >3, 


( 2 . 6 ) 


where jS'’* denotes the surface measure of the Euclidean sphere in e.g., |§^| = 27r 
and |§^| = dvr. The classical simple-layer potential Tl and the double-layer potential ^ are 
formally dehned by 

:= [ G{x,y)ij{y)dy, (kv){x) := j d^^y)G{x,y)v{y) dy, x^W^\dVL] 
Jdn Jan 

here, denotes the (outer) normal derivative with respect to the variable y. These 

pointwise dehned operators can be extended to bounded linear operators 

^ e and ^ e L{H^/\dn);H\U \dn)). (2.7) 

It is well-known that A'tSip = 0 = A^v in U \ dO, for all G and v G 

The simple-layer, double-layer, adjoint double-layer, and the hypersingular integral operator 

are dehned as follows: 

= + =-^ + 7 ^*^, and 2n =- 7 )“'^. (2.8) 

These linear operators are bounded linear operators for s G { — 1/2, 0,1/2} as follows: 


QJ G i7^/2+"(aG)), (2.9) 

1^ G (2.10) 

G L(i7-^/2+"(aG);i7-^/^+*(aG)), (2.11) 

2n G (2.12) 

The operators QJ, 5J, ^ will often be applied to functions in T^(r). Throughout the paper, 
we employ the convention that for // G T^(r) we implicitly extend by zero, e.g., 

TJ// means 01(770^2/), means QJ(£'o,r 2 /), and ^'ip means .^'(Eoy//). (2.13) 


An analogous extension is obviously used when .^, QU are applied to an n G 

Remark 2.3. Ellipticity o/QJ and QU is not used in our analysis of Theorem 3.1 and Corol¬ 
lary 3.2. In particular, there is no need to scale G to ensure diam(G) < 1 m 2D or to assume 
that T is connected. ■ 

2.3. Surface simplices and admissible triangulations. Fix the reference simplex 
Tref := {x G 0 < Xi,..., Xd-i, < 1}) which is the convex hull of the d vertices 

(0, Cl,..., e^-i} (“0-faces”). The convex hull of any j -|- 1 of these vertices is called a “j-face” 
of Tref. We call the (d — 2)-faces “facets” of Tref. 

We require the concept of regular, shape-regular triangulations 7/ of T. 
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Definition 2.4 (regular and shape-regular triangulations). A set Th of subsets ofT is called 
a regular triangulation of T if the following is true: 

(i) The elements T ^ Th are relatively open subsets of T and each T is the image of Tef 
under an element map 7 t : Tei T. The element map is assumed to be bijective 
and on Tei- 

(a) The elements cover UreTh T -T. 

(Hi) “no hanging nodes”: For each pair (T, T') & Th x Th, the intersection T fl T' is either 
empty or there are two j-faces f, f C dTef of Tei with j e {0,... ,d — 2} such that 
TnF = 7T(/)=7T'(/0- 

(iv) Parametrizations of common boundary parts of neighboring elements are compatible: If 
(/) T nT = 7t(/) = jT'if'), than 7 ^^ o 7 ^-/ : f' ^ f is an affine isomorphism. 

We call the images of vertices ofTe^ under the element maps nodes ofTh and collect them 
in the set Mh- The images of the {d — 2)-faces of Tei are called facets ofTh and collected in 
the set Th- For each T G Th, we set h(T) := diam(T) := \x — y\. 

A regular triangulation is called K-shape regular, if the element maps 'Jt satisfy the fol¬ 
lowing: 

(v) LetGT{x) := 7 r(a^)"''- 7 r(a^) G be the symmetric Gramian matrix ofjT- The 

triangulation is K-shape regular if for allT E Th the extremal eigenvalues Xmin{GT{x)) 
and XmaxiGrix)) of Gt{x) satisfy 

( h{TY X^UGt{x))\ ^ 

Sf VA™n(Gr(x)) ^ hiTf ) -''■ 


(vi) If d = 2, we reguire explicitly that the element sizes of neighboring elements are com¬ 
parable: 


h{T) < Kh{T) for all T, T with T n T' ^ 0. 


With each triangulation Th, we associate the local mesh size function h G L°°(r) which is 
defined elementwise by h\T := h{T) for all T G Th- We note that for a K-shape regular mesh 
we have 


max 

T&Th 


h{TY-^ 

|r| 


< 1, 


(2.14) 


where the implied constant depends solely on k. 

If r is the union of pieces of (d — l)-dimensional hyperplanes and the element maps 
are affine, then the Gramians are constants and the Definition 2.4 generalizes the classical 
concept of a shape-regular triangulation of F. In the non-affine case, the following example 
illustrates how triangulations as stipulated in Definition 2.4 can be created: 


Example 2.5. Let F C dQ be an open surface piece and assume F = 7 (F) for some reference 
configuration F C and some sufficiently smooth map 7 . Let Th = {T,... jTat} be a 

standard, regular, shape-regular triangulation ofV with affine element mapsfif,, i = 1,..., N. 
Then, the triangulation with elements T = 7 o 7 .^ (Tref) and element maps 7 o 7 ^ satisfies 
the hypotheses of Definition 2.4- 

This concept generalizes to surfaces consisting of several patches; it is worth emphasizing 
that in that case the patch parametrizations need to match at patch boundaries. ■ 

6 







Remark 2.6. In Definition 2.4 the conditions on the mesh are formulated so as to ensure 
that the spaces S^~^^{Th) below have good approximation properties. The conditions (Hi) and 
(iv) in Definition 2.4 could he relaxed if only good approximation properties of the spaces 
'P’^i'Th) are reguired. ■ 

For an element T & Th, 'WQ define the element patch uJh{T) by 

uJh{T) := (1J{T^ : T'G with T n ^ ^ 0})° . (2.15) 

The assumptions on the element maps of a K-shape regular triangulation imply that elements 
of a patch are comparable in size. Furthermore, the fact that F results from a Lipschitz 
dissection of dQ imposes certain topological restrictions on the patches: 

Lemma 2.7. Let Th be a regular, n-shape regular mesh. Then there is a constant C > 0 
that depends solely on k and the Lipschitz character of dQ such that the following holds: 

(i) h(T) < Ch(T') for any two elements T, T' with T n T' 7 ^ 0. 

(a) The number of elements in an element patch is bounded by C. 

(Hi) For any two elements T, T' in the element patch Uh{T") there is a seguence T = 
Tq, ... ,Tn = T' of elements Ti, i = 0,... ,n, in UhiT") such that two successive elements 
Ti, Tj+i share a common facet: Ti n Tj+i G IFh for i = 0 ,..., n — 1. 

Sketch of Proof. Statement (Hi): We show (hi) hrst for the node patch 

ujh{z) := {T : T eTh with 2 ; G T} j 

and some node z of T". This follows from the fact that F results from a Lipschitz dissection 
and considerations in using local charts. After a Euclidean change of coordinates, we 
may assume that dQ is (locally) a hypograph, i.e., there is a Lipschitz continuous function 
A : Rr(0) —)■ M such that the set {(a:,A(x)): x G 5^(0)} C dll. Without loss of generality, 
we assume the Euclidean coordinate change is such that z = (0, A(0)). One may also assume 
that A is dehned on (and Lipschitz continuous) so that the map A : —)■ given by 
(x, t) i-G- (x, A(x) + t) is bilipschitz. 

We distinguish the cases 2 ; G F and 2 ; G 9F. Let 2 ; be an interior point of F. Then, the pull¬ 
backs T := A“^(T), T C ijOh{z), are contained in the hyperplane x {0} and (identifying 
this hyperplane with completely cover a neighborhood of 0 G This together 

with (hi) of Dehnition 2.4 shows the claim. If z G (9F, then the fact that the elements are 
contained in F and that F results from a Lipschitz dissection implies that near 0 G 
the pull-backs T are all on one side of a Lipschitz graph in This together with (hi) 

of Dehnition 2.4 again implies the claim. Since uJh{T") is the union of the d node patches 
u!h{z) associated with the d nodes of T", this concludes the proof of (hi). 

Statement (ii): Consider the case of an interior point z G F. The assumption (hi) of 
Dehnition 2.4 and the fact that the map A is bilipschitz implies that the solid angles of the 
elements T at 0 are bounded away from zero by a constant that depends solely on k and A. 
This implies the claim for a node patch 0 Jh{z) and thus for Uh{T) with T G Th- 
Statement (i): For d = 2, this follows by dehnition. For d > 3 we hrst note that two elements 
sharing a facet f & Th have comparable size by (hi)—(v) of Dehnition 2.4. We conclude the 
proof with the aid of statements (hi) and (ii). □ 
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2.4. Admissible weight functions and discrete spaces. 

Definition 2.8 (a-admissible weight functions and polynomial degree distributions). A func¬ 
tion Wh G L°°{r) is (T-admissible with respect to Th if 

||'«^h||L°°(T) < o'Whix) almost everywhere on oJhiT). 

A a-admissible function G L°°(r) is called a a-admissible polynomial degree distribution 
with respect to Th, if QhiT) := gtlr G No for all T G Th- 

We write 

T^{Th) := G L^(r) : W eTh o 77 . is a polynomial of degree < qh{T)}, (2.16) 

for the space of (discontinuous) piecewise polynomials of local degree qh{T). Moreover, we 
introduce spaces of continuous piecewise polynomials of local degree qh{T) -|- 1 by 

S^+^Th) := P^^^Th) n H\T), (2.17) 

S^+^Th) := S^+^Th) n H\T). (2.18) 

We note the inclusions V^{ThlC L‘^{T) C H-P^T), S^+\Th) C #^(r) C HP^{T), and 
S^+\Th) C H\T), as well as = S^+^T) in case of h = dn. 

For q G No, the use of non-boldface superscripts in V^iTh), S^^^iTh), and S^^^^iTh) indi¬ 
cates that a constant polynomial degree is employed. 

3. Main result and applications 

3.1. Inverse estimates. The following Theorem 3.1 is the main result of this work. 

Theorem 3.1. Let Th be a regular, n-shape regular triangulation ofV and let Wh G L°°(r) 
be a a-admissible weight function with respect to Th- Then, it holds 

||'«ihVrTJ'0||i2(r) -I- \\wh^''f>\\L‘2(r) < C'inv(||'«ih/h^'^^||L°°(r)||'0llH-i/2(r) + (3.1) 

||'W^ftVr.fiT||L2(r) -|- \\whiL!0v\\L2^-p) < C\m{\\'0Jh/h^^‘^\\L<=°{T)\\v\\Hi/2^Y) + ll'*^/t^rT||L2('r)), (3.2) 

for all functions fj G T^(r) and all v G 77^ (T). The constant Cinv > 0 depends only on dQ, 
T, the K-shape regularity ofTh, and a. 

In the following Corollary 3.2, we apply the estimates (3.1)-(3.2) of Theorem 3.1 to discrete 
functions T/i G V^iTh) and I 4 G S'^^^iTh)- We mention that the restriction to <7 G {2,3} 
in Corollary 3.2 is due to the fact that the underlying reference [KMR14] restricts to this 
setting. 

Corollary 3.2. Let Th be a regular, K-shape regular triangulation of T. Suppose that d G 
{2,3} and that qh is a a-admissible polynomial degree distribution with respect to Th- Then, 
there exists a constant Cinv > 0 such that the following estimates hold: 

\\hP\qh + 1)-' Vr5J^/.||L2(r) + \\hP\qh + 1)“' < anv||^/^b-i/2(r), (3.3) 

\\hP\qh + 1)-'VrJ^WlU^CD + + 1)"'2IJW|U2(r) < anv||14bi/2(r), (3.4) 


for all discrete functions G V^{Th) and Vh G The constant Cjn^ > 0 depends 

only on dQ, F, the K-shape regularity of Th, and the a-admissibility of q^, but is otherwise 
independent of the polynomial degrees and the mesh Th- 

Proof. The starting point are the following two inverse estimates 

\\h^/\qh + 1)-' ^hh^r) < ll^h||/^-i/ 2 (r) for all T, G (3.5) 

+ l)-iVrfo.||L^(r) < imsmiD fo^ ah G (3.6) 

where the hidden constants depend solely on dQ, F, the K-shape regularity of Th, and the 
cr-admissibility of qh- The bound (3.5) is essentially taken from [GeoOS, Thm. 3.9]. However, 
since the non-trivial interpolation argument is not worked out in [GeoOS, Thm. 3.9] and 
since [GeoOS, Thm. 3.9] is not concerned with open surfaces F, we present the details in 
Lemma A.l. We remark that its proof employs the characterization of fractional Sobolev 
norms in terms of the Aronstein-Slobodeckii norm. The bound (3.6) follows also from poly¬ 
nomial inverse estimates and an interpolation argument for spaces of piecewise polynomials, 
which is non-trivial—see [KMR14] for details. We also refer to [AFF^14, Proposition 5] for 
the h-version of (3.6), in which the dependence on the polynomial degree qh is left unspecihed. 

We dehne a weight function by Wh ■= h^^‘^{qh + 1)”^- Note that \\wh/h^^^\\L°°{r) < 1 
and that Wh is r-admissible, where r depends only on k and a. The combination of (3.5) 
with (3.1) leads to (3.3). The bound (3.6) in conjunction with (3.2) yields (3.4). □ 


3.2. Application to efficiency of residual error estimation. 

3.2.1. Weakly singular integral equations. The next corollary proves that the estimate (3.1) 
provides stability of 05 and in weighted norms for subspaces (1 — P/i)L^(F) C T^(F), where 
Ph is some projection operator. Note that the following corollary is in particular applicable 
to the Galerkin projection onto V^{Th)- 

Corollary 3.3. Let Th be a regular, K-shape regular triangulation ofV. Let Xh be a closed 
subspace of with V^{Th) P Xh C T^(F). Let H/^ : L^{T) —)■ Xh be the Lf-orthogonal 

projection onto Xh and P/i : iF“^/^(F) —)■ Xh C iF“^A(F) denote an arbitrary 
stable projection onto Xh. Then, there is a constant Cinv > 0 depending only on the K-shape 
regularity ofTh, the stability constant ofFh, on dQ as well as F such that for all (f G T^(F) 
and Ph G 

]|/,i/2vr55(i - Ph)ct>\\L^(v) + \\h^/^^{l - Ph^mv) < - Ph)(l>\W(T). (3.7) 

Proof. For arbitrary w G H^P{dkl) we get by transformation to the reference element and 
standard approximation results that ]|(1 — n/i)tc]|^ 2 ( 7 -) ^ ^(7")l|wl|^i/2(7^), where we employ 
the Aronstein-Slobodeckii norm in the definition of ]| • \\hi/^{t)- Hence, by summation over 
all T G Th, using the Aronstein-Slobodeckii characterization of ]| • ]|Hi/ 2 (an), and then the 
characterization (2.1) of the norm ]| ■ ]|iifi/ 2 ('r), we arrive at 

]|/i-V 2 (i _ n;j)t(;]|^ 2 (r) < ]|M;]|j;^i/ 2 (r) for all w G 
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Orthogonality of 11^ and a duality argument then shows (see [CP06, Theorem 4.1] for the 
analogous proof on polygonal boundaries.) 

11(1 - n,)0||;^_v2(r) < ||/^'/"(l - n,)0|U2(r) for all 0 G L^T). (3.8) 

Combining this estimate with the inverse estimate (3.1) for tjj = (1 — T\.h)4> Wh = 
we get 

||/,i/2VrQJ(l - n,)0|U.(r) + - n,)0|U.(r) < for all 0 G ^'(r). 

For an if“^/^(r)-stable projection P/j, we note that the projection property of P/j implies 
(1 — P/i)(l — II/i) = (1 — P/i). This and elementwise stability of 11^ imply, for all 0 G T^(r), 

11(1 - P»)<^|ls-.„,r, < 11(1 - ||A'''^(1 - n»).(.||y,r, < 

Finally, we use the projection property (1 — P/^)^ = (1 — P/^) and argue as for II/i to obtain 

||/,i/2VrQj(i - P,)0||i.(r) + - P/^)0 ||l2(p) < ||h'/2(l - P;,)0|U2(r) for all 0 G ^^(F). 

This concludes the proof. □ 

One immediate consequence of Corollary 3.3 is the efficiency of the weighted residual error 
estimator rjh from [Car97, CMSOlj: Suppose that 03 is i7“^/^(F)-elliptic (in the case d = 2, 
this can be enforced, for example, by the scaling requirement diam(O) < 1). For / G ih^(F), 
let 0 G be the unique solution of the weakly singular integral equation 030 = /. 

Let Xh C T^(F) be a discrete space which contains at least the piecewise constants, i.e., 
V^iTh) ^ Xh, and let G Xh be the unique Galerkin approximation of 0 in Xh, i.e., 

(03(0 - $;,), T;,)r = 0 for all ^>h e Xh- (3.9) 

Under these assumptions (and, strictly speaking, for polyhedral F), [CMSOl] proves the 
reliability estimate 

C,;,‘ 110 - < >0.® ■= ll'!‘'"Vr(/ - 9J4»)lli>(r). (3.10) 

The constant Crei > 0 depends only on F and the K-shape regularity of Th- The following 
corollary provides the converse efficiency estimate with respect to some slightly stronger 
weighted L^-norm. We note that the additional assumption 0 = 03“^/ G T^(F) is in partic¬ 
ular satished for F = dVL. 

Corollary 3.4 (Efficiency of 7]h,‘^ for weakly singular integral equations). Let Th be a regular, 
K-shape regular triangulation o/F. Assume 0 = G T^(F) and let Xh C iL“^/^(F) be 

a closed subspace with V^{Th) ^ Xh C T^(F). Let G Xh be given by (3.9). Then the 

weighted residual error estimator from (3.10) satisfies 

Vh,<n < Ces ||^^'^^(0 — ®/i)||L2(r), (3-11) 

where = Cinv > 0 zs the constant from Corollary 3.3. 

Proof. With the Galerkin proiection Pi, ; —)■ Xh and $!, = P(j0, Corollary 3.3 yields 

Vh,^ = - <^>h)llLhr) < " ®/^)IU2(^). □ 
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Remark 3.5 (Stronger efficiency of 2D BEM). While the efficiency estimate (3.11) involves 
a slightly stronger norm on the right-hand side, particular situations (as, e.g., the 2D direct 
BEM formulation of the Dirichlet problem [AFF'*’13b]j permit to bound ||/i^/^(0 — 
by up to higher-order terms. In [AFF’'‘13b], this is achieved by decomposing 

(j) in a singular part associated with the vertices of D and a regular part; the higher-order 
terms depend only on the regular part of (j). ■ 

3.2.2. Hypersingular integral equations. Results similar to Corollary 3.3 also hold for the 
double-layer integral operator ^ and the hypersingular integral operator OR. Here, particu¬ 
larly interesting choices for the projection P/^ are Scott-Zhang type projections onto 
see [SZ90] as well as the adaptation to BEM in [AFF’''14, Section 3.2]. 

Corollary 3.6. Let Th be a regular, n-shape regular triangulation ofT. Let Xh C be 

a closed subspace with S^{Th) C Xh C H^{T). Let P/j : —>■ Xh be an H^^^{T)-stable 

projection onto Xh. Then, for all v G R^(r), 

\\h^/^VrIi{l - Fh)v\\L2ir) + ||h'/'2n(l - Fh)<j\\LHT) < - Ph)w||L 2 (r). (3.12) 

The constant Cinv > 0 depends only on the K-shape regularity of Th, the stability constant of 
P/i, and r. 

Proof. Arguing along the lines of the proof of Corollary 3.3, we hrst consider the Scott-Zhang 
projection Jh : —)■ S^{Th) onto S^{Th). According to [AFF’'‘14, Lemma 7] (strictly 

speaking, this result is formulated for polygonal boundaries only, but the proof transfers 
with minor changes to the present case), it holds 

11(1 - 4)u;bv2(r) < min \\h^/^V{w - Wh)\\L^ir) < - 4)u;|U=(r) (3.13) 

Wh&S\Th) 

for all w G H^{T). The hidden constant depends only on T and the K-shape regularity of Th- 
Combining this with the inverse estimate (3.2) for u = (1 — Jh)w, we arrive at 

\\h}/^Xrm-Jh)w\W;r) + \\h^'^W{l-Jh)w\\L 2 ;r) < \\h^l^Vr{l-Jh)w\W;r) for all w G H\T). 

The same arguments apply for any iL^/^-stable projection P/j ; —)■ Xh, but addition¬ 
ally employ its stability. □ 

As in Section 3.2.1, an immediate consequence of Corollary 3.6 is the efficiency of the 
weighted residual error estimator rjh from [Car97, CMPS04] for the hypersingular integral 
equation: Suppose that 77^/^ (T) does not contain any characteristic function Xuj for cu C T 
(this is in particular satisfied if dD is connected and T ^ dD). Then, 211 : 77^/^ (T) —)■ 
77 - 1 / 2 (T) jg elliptic isomorphism. For / G L^(r), let u G 77^/^ (T) be the unique solution 
of the hypersingular integral equation ORu = /. Let Xh C 77^ (T) be a discrete space which 
contains at least the piecewise affines, i.e., S^{Th) C Xh- In addition, let Uh G Xh be the 
unique Galerkin approximation of u in Xh, i.e., 

(QR(u - Uh), Vh)r = 0 for all Vh G Xh- (3.14) 

Under these assumptions (and, strictly speaking, for polyhedral T), [CMPS04] proves the 
reliability estimate 

C,;,‘ ll« - t'j.k./.iD < -■= \\h''^u - <muo\\LHT). 
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(3.15) 


The constant C^ei > 0 depends only on T and the K-shape regularity of Th- The following 
corollary provides the converse efficiency estimate with respect to some slightly stronger 
weighted TT^-seminorm. 

Corollary 3.7 (Efficiency of ri^ for hypersingular integral equations). Let u = 211“^/ G 
H^{T). Let Xh C be elosed with S^{Th) ^ Xh C H^{T). Let Uh G Xh he given by 

(3.14). Then the weighted residual error estimator from (3.15) satisfies 

Vh,w < Ces ||h^'^^Vr(M — Uh)\\L'2{r), (3.16) 

where = Cinv > 0 zs the constant from Corollary 3.6. 

Proof. With the Galerkin projection —)■ X^ and Uh = FhU, Corollary 3.6 yields 

Vh,w = ||hV22n(n - Uh)\\mr) < \\h^/^Vriu - Uh)\\mr). □ 

Remark 3.8. //T = dfl is connected, the kernel o/QU is the space f constant functions on 
T. Therefore, 211 : —)■ is an elliptic isomorphism, where H^{dkl) := 

G : {v, l)^^ = O} for |s| < 1. Recall that 211 : Hl{dPl) —)■ Hl~^{dPl) is 

an isomorphism for all 0 < s < 1. For f G H^{dQ), the solution u := 211“^/ thus has 
additional regularity u G Hl{dVL), and Corollary 3.7 holds accordingly. 

3.2.3. Remarks on the extension to hp-BEM. The above efficiency statements are formulated 
for the h-version BEM. They do generalize to the hp-version. 

Since the corresponding reliability estimates have only been formulated for closed surfaces 
T = dQ and affine element maps in [KM14], we restrict the following result to this setting: 

Corollary 3.9. Let d G {2,3}, T = dQ, and Th a regular, K-shape regular triangulation of 
T. Assume that the element maps are affine. Let Qh he a a-admissible polynomial degree 
distribution. Then there exists C > 0 depending only on dQ, the K-shape regularity of Th, 
and a such that the following holds: 

(i) Let 2J : —)■ be an isomorphism. Let 0 = 2J“^/ for some f G 

Set Xh := V^iTh) and let ^hp G Xh be the Calerkin solution given by (3.9). 


Then: 

C ^110 — 4’/ip||/f-i/2(aQ) < Php.Qj := ||(^/(1 + ?/x))^'^^Vr(/ — 2J$/ip)||p,2(gQ), (3.17) 

Vhp,<n < C\\{h/{1 + qh)y^‘^{4> — ^hp)\\L'^{dn)- (3.18) 

(a) LetT be connected andu = W~^f for some f G H^{dLl). Set Xh := S'^~^^{Th)CHl{d^l) 
and let Uhp G Xh be the Calerkin solution given by (3.14). Then: 

C ^||m — <3hp,w'-= \\{h/{1 + qh)Y^‘^{f — ^Uhp)\\L'^{dn)-, (3.19) 

Vhp,w ^ C \\{h/{l 3-qh)y^‘^Xr{u — Uhp)\\L^{dn) (3.20) 

+ II(^/(l + — Uhp)\\L^{dn) . 

Proof. The reliability bounds (3.17), (3.19) are taken from [KM14, Cor. 3.9, Cor. 3.12]. 


For the proof of (3.18), we let Ilhp be the L^(cir2)-projection and be the Calerkin pro¬ 
jection. We note that the analogue of (3.8) is 

||(l-P,p)0||^-, 

/2(ao) ~ 11(1 - Thp)t/j\\H-i/2^9n) < ||(V(1 + V0 G L'^{dn) 

(3.21) 
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(cf. [KM14, Thin. 3.8]). Hence, proceeding as in the proof of Corollary 3.3 with Wh = 
(V(l + qh)Y^^ we get 

— P/ip)0||L2(aQ) ^ 11(1 + qh) ^'^^||l°°( 90)||(1 — P/ip)^0||H-i/2(ao) + ll'W^fe(l “ Ph)0||L2(ao) 

( 3 . 21 ) 

~ 11(1 + qh) ^'^^||L°°(aO)||'«^/i(l — P/ip)0||L2(aQ) + ||w/i(l — P/ip)0||L2(aQ). 

The proof of (3.20) proceeds along similar lines. The key is the analog of (3.13). Combining 
[KM14, Lem. 3.7] and the proof of [KM14, Lem. 3.10] produces an approximation operator 
J'p ; H\dQ) ^ S^^\Th) with 

11(1 “ Jhp)'^\\H^/^{dQ.) ~ ll(^/(l + qh)Y^‘^'^rv\\L'^{an) + ||(^/(1 + <?h))^'^^'«^||L2(ao)- 

Finally, an operator J^p : Hl{dfl) —?■ X^p is then obtained by setting J^pV := J'f^pV — 
J'hpV, where the overbar denotes the average over dQ. It is easy to see that Jhp has the 
same approximation properties as on the space Hl{dQ). Proceeding as in the proof of 
Corollary 3.3 or 3.6 hnishes the proof. □ 

4. Far-field and near-field estimates for the simple-layer potential 

The proof of Theorem 3.1 is based on decomposing the pertinent potentials into “far-held” 
and “near-held” contributions. In the present section, we analyze the decomposition for 
the simple-layer potential and provide inverse estimates for both components. Section 4.2 is 
concerned with inverse estimates for the near-held parts, which essentially follow from scaling 
arguments, whereas Section 4.3 deals with the far-held part. Throughout the section, we let 

"0 G L^(F) and assume that ^|J is extended by zero to Oil \ F, (4-1) 

i.e., we identify ^|J with 

4.1. Decomposition into near-held and far-held. For a parameter 5 > 0, we dehne 
for each element T ETh the neighborhood Ut of T by 

T C Ut ■= [^ B2Sh{T){x)- (4.2) 

x£T 

Since dQ is Lipschitz and F stems from a Lipschitz dissection and by K-shape regularity of 
%, we can hx the parameter <5 > 0 and hnd M G N (both 6 and M are independent of Th) 
such that the following two conditions are satished: 

(a) F n Ut is contained in the patch uJh{T) of T (see (2.15) for the dehnition), i.e., 

TnUT^UhiT). (4.3) 

(b) The covering F C UtgT,, is locally hnite with a uniform bound, i.e., 

#{17t : T G T; and a: G Ct} < M for all a: G (4.4) 

Finally, we hx a bounded domain 17 C such that 

UtCU for all T G 77. 
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(4.5) 



It will be important that U is chosen independently of Th- To deal with the non-locality of 
the integral operators, we define for T E Th the near-held and the far-held 
simple-layer potential := fO'ip by 

:= ^(V’Xrnc/r) and ■= ^(TXtxut)^ (4-6) 

where Xi^ denotes the characteristic function of the set u: C We have the obvious identity 

forallTeTft. (4.7) 

In our analysis, we will treat and separately, starting with the simpler case of 


4.2. Inverse estimates for the near-held part The near-held parts of a potential 

can be treated with local arguments and the stability properties of the associated boundary 
integral operators. 


Lemma 4.1. There exists a constant Cnear > 0 depending only on dQ, F, and the n-shape 
regularity ofTh such that for arbitrary T & Th and G V^iTh) with supp C ojh{T) it 
holds 

iivS'i-riu^y/T) < 


Proof. We hx an element T ^ Th- We recall that is piecewise constant and compute for 

X & Q 

(V^T^)(a:) = / ^xG{x, y) dy for all a: e \ T. 

T'euJhiT) 


The number of elements T' in the patch Uh{T) is bounded in terms of the shape regularity 
constant k (cf. Lemma 2.7). With some constant that depends only on k, we bound 


|(VQJTD(x)P 


< 


E 

T'&uJhiT) 




h\T'\ 


\T ^G{x,y)\dy 


IT’ 


(4,8) 


Next, we show for elements T' C Uh{T) 

j (j \V.G{x,y)\dy)" dx<h{TY. (4.9) 

J JJ'jp J T' 

This follows from a local Lipschitz parametrization of dVt. We assume that (after possibly a 
Euclidean change of coordinates) that {(x', A(x')): x' G i?2r(0)} is a part of dVL that contains 
Uh{T). The function A is Lipschitz continuous, and we remark in passing that by [Ste70, 
Thm. 3, Sect. VI] we may assume that A : —)■ M is Lipschitz continuous. (If such a local 

consideration is not possible, then, since the number of local charts is finite by definition of 
bounded Lipschitz domains, we must have diam(a;/i(T)) = 0(1) so that (4.9) is trivially true.) 
We may also assume that Ut C {(x', A(x') -1- f) | x' G i?2r(0), t G M}. The key observation is 
that the mapping A : ^ given by {x',t) ^ (x', A(x') -|- t) is bilipschitz. We conclude 

that A~^Ut ='■ Ut and A~^T' =: T' C Br(0) x {0} satisfy for some xq and some c > 0, which 
depends solely on the bilipschitz mapping A, 

Ut C Bch{T){.Xo) X [-ch(T),ch(T)], f' C Bch{T){.Xo) x {0}. 
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Finally, using VxG{x, y) ~ \x — y\ the definition of the surface integral, and the change 

of variables formula for bilipschitz mappings from [EG92, Sec. 3.3.3], we get 


' x&Jt 


/ \VxG{x,y)\ dy\dxc:^ l 
lyeT' J JxeUr \JyeT' 


< 


X — y\ dy'^ dx 


r rch{T) 

?GSch(T)(a:o) Jt=-ch(T) \-'r;GBch(T)Go) 

d I I if „|2 I .2\-('^-1)/2 


dy ) dtd^ 


^ h{Tr I / / (I? - >)P + G 

J^€Bc{xo) Jt=—c \JriGBc{xo) 

~ h{Ty, 


dy ) dtd^ 


where the last estimate follows by a direct estimation of the integrals, which is independent 
of h{T). We have thus shown (4.9). Inserting (4.9) in (4.8) gives 



\{vm<!/l)(x)\^dx < Y. K\T,\^h(Tr~ 

T'&u>h{T) 




□ 


Proposition 4.2 (Near-field bound for QJ). Let Wh be a a-admissible weight function. There 
exists a constant G^ear > 0 depending only on dVL, F, the n-shape regularity of Th, and a, 
such that the near-field part satisfies G H^{U) and G H^iX) together with 

Y l|tt').Vr7;”‘«“lli.(T) + Z ll«''./'‘‘''"lli-(T)l|V«S’,¥111.(1,,) < C.e„ (4-10) 

Ten Ten 

Proof. The stability (2.9) of QJ ; proved in [Ver84] gives, for each T E Th, 

IIVr7o“*'*^»,TlU2(r) < \\^{i’XUTnr)\\m(dn) ^ ||V’X;7Tnr||L2(aQ) = ||■^/’||L2(c/J,n^)• 

Summing the last estimate over all T G 7/i and using (4.3)-(4.4), and a-admissibility of Wh, 
we arrive at 

ll'^/i^r7cr*'*^»(TII l2(t) ^ ll'^/illL°°(T)lll('llL2(;7j,nr) — ll^/it/'llL2(r)5 
Ten Ten 

where all estimates depend only on the K-shape regularity of Th and the admissibility constant 
a. This bounds the first term on the left-hand side of (4.10). To bound the second term, 
let fl/j denote the L^(F)-orthogonal projection onto V^{Th)- We decompose the near-field 
as = ^{UhifjXrnUT)) + ^((1 “ ^h)i’XrnUT)- The condition supp(V’Xrn(7T) ^ (^h{T) 

implies supp {Th{fJXrnUT)) ^ ^h{T) and therefore, taking = lihii’XvnUT) Lemma 4.1, 
we conclude 

Ten 

~ ll'^/i/^^^^llL°°(T)ll^^^^n/*(V'Xrni7T)llL2(^^(7,)) (4T1) 

Ten 

^ Wwh'ifWl^iry 
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where we used the local L^-stability of II/i in the last estimate. Recalling the stability 
QJ : —)■ H^{U) of (2.7), the equality (2.4), and the approximation property (3.8) 

of fl/i, we get 

T&Th 

( 2 . 7 ) 

~ “ n/i)'0Xrni7Tll//-i/2(r) (4-12) 

Ten 

^ ll^/^/^^^^lli-(r)ll^^^^(V'Xrnc/T)lli2(r) - lkhV'lli2(r). 

Ten 

Combining (4.11)-(4.12), we bound J^Ten ll'^^/^^'^^lli°°(T)ll^'^»,TllL 2 (i 7 y) conclude the 
estimate in (4.10). □ 

4.3. Estimates for the far-held part The following lemma is taken from [FKMP13]. 

For the convenience of the reader and since the same argument underlies the proof of the 
analogous lemma for the double-layer potential (Lemma 5.3 below), we recall its proof here. 

Lemma 4.3 (Caccioppoli inequality for With 12®^* = \ 12, the function from 

(4.6) satisfies G C°°(12), g C°°(12®'^*), and G C°°{Ut)- Moreover, 

there exists a constant Cca.cc > 0 depending only on dVL, F, and the n-shape regularity of Th 
such that Hessian matrix satisfies 

11-^ ^®,rlU^(B 5 j,(T)(a:)) — ^cacc \\'^'^T},T^\L^{B 2 Sh{T)i^)) all X ^ T E Th- (4-13) 

Proof. The statements G C°°(12) and g C'°°(12®'^*) are taken from [SSll, 

Theorem 3.1.1], and we therefore focus on ^ C°°{Ut) and the estimate (4.13). 

According to [SSll, Proposition 3.1.7], [SSll, Theorem 3.1.16], and [SSll, Theorem 3.3.1], 
the function G := |n : —)■ M : v\k G H^{K) for all K <Z compact} 

solves the transmission problem 

= 0 in 12 U 12®^*, 

[n|y = 0 in 27^/2 (912), (4.14) 

bi^nr] = -bXv\UT in 27"^/^ (912). 

In particular, (4.14) states that the jump of as well as the jump of the normal de¬ 
rivative vanish on 912 fl Ut- This implies that u<xs,t is harmonic in Ut by the following 
classical argument: First, we observe that is distributionally harmonic in Ut-, since a 
two-fold integration by parts that uses these jump conditions shows for v G CT{Ut) that 
(n® T) = 0. Weyl’s lemma (see, e.g., [Mor08, Theorem 2.3.1]) then implies that 

is therefore strongly harmonic and G C°°{Ut)- 

The Caccioppoli inequality (4.13) now expresses interior regularity for elliptic problems. 
Indeed, for each u G H^{Br+e) such that u G and Au = / on Bj.+s with balls 

Br C Br+e with radii 0 < r < r -f £ and some / G L‘^{Br+e), [Mor08, Lemma 5.7.1] shows 

||-D^'u]|l 2 ( 5 ^) < (^\\f\\L^[Br+e) + “ II V-u]|i, 2 ( 5 ^_^^) ^||'u||L2(n,+,)y ( 4 - 15 ) 
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The hidden constant depends solely on the spatial dimension and is independent of r,£ > 0, 
and u, f. We apply (4.15) with r = 6h(T) = e, f = 0, and u = — ct, where ct = 

\B 2 Sh(T){x)\ ^B 2 Sh(T){x) ^y- additional Poincare inequality finally leads to (4.13). Note 

that 5 and hence (Peace depend only on dVL, P, and the K-shape regularity of Th- □ 

The non-local character of the operator QJ is represented by the far-held part. Lemma 4.3 
allows us to show a local inverse estimate for the far-held part of the simple-layer operator; 

Lemma 4.4 (Local far-held bound for QJ). For all T G Th, it holds 

||hl/2Vr7i“*<TlU2(T) < l|/i'/V<TlU2(T) < Qar ||. (4.16) 

The constant (Pfar > 0 depends only on P, and the n-shape regularity constant ofT- 

Proof. By Lemma 4.3 we have ^ C°°{Ut). The hrst estimate in (4.16) follows from the 
fact that, for smooth functions, the surface gradient Vr(-) is the orthogonal projection of 
the gradient V(-) onto the tangent plane, i.e., Vr 7 o“*M(x) = Vu{x) — (Vm(x) ■ ^'(x), 

see [Ver84]. 

The second estimate in (4.16) is proved with a trace inequality and the Caccioppoli in¬ 
equality (4.13) in the following way. We hx an element T G Th- 

1. step: We provide a trace inequality. Let B = Br{x) be a ball with center x G T C dQ and 
radius r > 0. Let B' = i? 3 r/ 2 (x) and B" := B^r/A{x). We dehne a smooth cut-oh function 
Xb £ (P“(M‘^) with supp Xb F B' and = 1 on i? by 

Xb ■= Xb" * Pr I a, 

where Ps{x) is a standard molliher of the form Ps{x) = e~'^pi{x/e) for a hxed pi G (P“(M'^) 
with Pi > 0, supp Pi C i?i(0) and pi{x) dx = 1. We note that for a (P > 0 depending 
solely on the choice of pi, we have 

llVxBllL^(Md) < (Pr"h 

With this cut-off function in hand, we estimate for sufficiently regular functions v and the 
standard multiplicative trace inequality for dQ 

\\v\\h(Bndn) < 11X5^11x2(911) < \\xbv\\12(^q) + \\xbv\\l^(q)\\V{ xBv)\\L^(n) 

^ + MmB')\\^v\\L2{B')- (4.17) 

2. step: The set T ;= ^^Bsh{T)/ 2 {x) | x G t| is a closed cover of T with sup^gj- diam(i?) < oo, 

and T is the set of their midpoints. According to Besicovitch’s covering theorem, cf. [EG92, 
Sect. 1.5.2], there is a constant AP, which depends only on the spatial dimension d, as 
well as countable subsets Qj F T, j = 1,... ,Nd, the elements of every Qj being pairwise 
disjoint, such that T C IJ^i Uxeg Let Qj be the set of balls obtained by doubling the 
radius of the balls of Qj, i.e., Qj := [Bsh(T){x) \ Bsh(T)/ 2 {x) G Gj}- As the elements of Qj are 
pairwise disjoint and all balls have the same radius 6h{T)/2, there is a constant N^, also 
depending only on the spatial dimension d, such that at most elements of Gj overlap. If 

we write B := Bsh{T)/ 2 {x), B' := B 3 / 4 ^sh{T){x), and B := Bsh(T){x), the multiplicative trace 
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inequality (4.17) and the Caccioppoli inequality (4.13) show 


(4.17) 1 (4.iaj 1 


3. step: We use the last estimate to get 

Nd 


Nd 


. far 11 2 


< 


j —1 


T\\L^{BnT) rs^ 

•?'=1 Beg 


H E iiv“§; 


far ||2 <- ^dNd ||^_far ||2 

,rll ^2rR^ ~ TTT^ II “»,tIIl2(;7^). 


Il 2 (b) ^ 


This concludes the proof of (4.16). 


□ 


Summation of the elementwise estimates of Lemma 4.4 yields the following result: 

Proposition 4.5 (Far-held bound for TT). There is a constant Cfar > 0 depending only on 
dVL, r, the K-shape regularity of Th, and the a-admissibility of the weight function Wh such 
that 

||«.ftVr7;“«g'Tlli.,T) < E h’l>^<T\\h(T) 

Ten Ten 

< Cfar (lkhA^^^||i-(r)||V'll|-i/2(r) + W^hHl^r)) ■ 

Proof. We use the local far-held bound (4.16) of Lemma 4.4 and to see 

(4.16) 

ll'^/i^r7(r*'*^®(TllL2(r) — II'*^/i/^^^^IIl°°(t) ll^'^®(rllL2(;7T) 

Ten Ten Ten 

~ ll'^fe/^^'^^lli°°(T)Il^^'l(’lli2(^.^) + |||2({;y). (4-18) 

Ten Ten 

The hrst term on the right-hand side in (4.18) is estimated by stability of QJ, the hnite 
overlap property (4.4), and (2.4) 

II'^/i/^^'^^IIl°°(t) ll^^l/^llL2(f/y) ^ ||'M^/i/^^'^^|lL°°(r)ll^^l/'lli2(f/) ^ ||'M^/i/^^'^^llL°°(r)llV'll|^-i/2(p) 

Ten 

The second term in (4.18) is bounded with the near-held bound (4.10). □ 


5. Far-field and near-field estimates for the double-layer potential 

Section 4 studied far-held and near-held estimates for the simple-layer potential. Correspond¬ 
ing results for the double-layer potential are derived in the present section. Throughout this 
section, let 

V e H\T) C H\dn). 

In particular, v G i7^/^(F) with ||u||/)^i/ 2 (r) = ||T||jpi/ 2 (gQ). Since does not allow jumps, 

the splitting into near-held and far-held contribution of the double-layer potential u^-.= 
cannot be done by characteristic functions, but requires smoother cut-oh functions and 
greater technical care. 
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5.1. Decomposition into near-field and far-field. We use the notation introduced in 
Section 4.1 concerning the neighborhoods Ut- In order to dehne the near-held and far-held 
parts for the double-layer potential, we need appropriate cut-oh functions: For each T G Th, 
we dehne 7]t G with the aid of the standard molliher that was already used in 

the proof of Lemma 4.4: 


Vt Ps/4h{T), 

Ux :— ^xeTBs/2h{T), 


(5.1) 

This function satishes: 




suppvT^Ux, Vt\u;^ 

III 

8 

ft. 

lA 

||Vr7T||L->(R4) < 

(6.2) 


where the implied constant depends on the K-shape regularity of the triangulation through 
the parameter 5. We note that the assumptions on Ut imply (suppr/^) fl F C uJh{T). 

The following lemma may be viewed as an extension of [DS80, Thm. 7.1] to the case of 
curved elements. 

Lemma 5.1 (Poincare-Friedrichs inequality on patches). Let v G For each T E Th, 

there is a constant vt EM. such that {y — Vt)pt £ {v — nT)(l — Vt) £ H^{dQ), and 

||n - VT\\L^{u,h{T)) < Ci\\hVrv\\L2(uih{T)), (5-3) 

||(n - VT)r]T\\H^/ 2 yn) < Ci\\h^''^Vrv\\L' 2 {ujHiT)), (5.4) 

||(n — VT)riT\\H^{dQ) < C'lllVr'y||L 2 (i.j;i(T)). (5-5) 

The constant vt satisfies vt = 0 if duh{T) fl dT contains a facet of the triangulation. The 
constant Ci > 0 depends only on dQ and the K-shape regularity constant ofTh- 

Proof. It is clear that {v — nT)(l — Vt) £ H^{dQ), since v E H^{r) and vt is smooth. The 
remaining statements require more care. 

1. step: For v G H^{r) and a facet f E Th oi the triangulation Th (recall that facets are 
images of {d — 2 )-faces of Tei under the element map) denote by if{v) the average of v 
on /. As if{l) = 1, we can use the Deny-Lions lemma on the reference element, and the 
assumptions on the element maps then imply 

W'V - £f{v)\\LyT) ^ h{T)\\Vrv\\LyT) if / is a facet of T, (5.6) 

^ if /i, /2 are two facets of T. (5.7) 

2. step: Fix an element T E Th- 

• If VT\dv = 0, then select an arbitrary facet fx of the element patch ujh{T). 

• If Vrldr ^ 0, then we claim that there exists a facet / of 0 Jh{T) with / C 5F. To see 
this, let xq E dr with Vt{xo) 7 ^ 0. By continuity of vt and since 9F is covered by 
facets of the triangulation, we may assume that xq is in the interior of a boundary 
facet fx- This facet belongs to a unique element T/ of the triangulation; by continuity 
of Vt, we may assume suppr^T fl T/ 7 ^ 0. Since (suppr^y) H F C Uh{T), we conclude 
Tf C Uh(T) and thus the boundary facet fx is a facet of UhiT). 

Set vx '■= ifj,{v). An immediate consequence of n G H^{r) is that uy = 0 if vt does not 
vanish on dT. Since vt is smooth, we conclude (n — vx)vt £ In fact, viewed as a 

function on dQ, we have 

supp((n - vx)vt) C uJhiT). 
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(5.8) 






3. step: The bounds (5.6), (5.7) in conjunction with Lemma 2.7 imply 

||n - VT\\L^{ujh{T)) ^ \\hVrv\\L^[ujh[T)), (5.9) 

II Vr(u - VT)\\L2{ujf^{T}) = II Vrv||L2(^^(r)), (5.10) 

where (5.9) is already the claim (5.3). The product rule, (5.2), (5.8), the estimate (5.9), the 
trivial bound h{T) < < 1 yield 

||Vr((v-t’r)hT)||L 2 ( 90 ) < Wiv - VT)'^rVT\\L^(ojhiT)) + llhTVr(v - t’r)||L2(a;ft(T)) < ||Vrr’||L2(<.j^(T)), 

which proves (5.5). It remains to verify (5.4). To that end, we recall the interpolation 
inequality ||M||^i/ 2 (gf^) < ||M||L2(an)||M||Hi(ao) for all u G H^{dn). Since |l(n - nT)? 7 T||i, 2 (an) < 

||t^ -'yT||L 2 (^^(r)), we get 

||(n - VT)VT\\H^/2(^9n) < IK^ - '>Jt)Vt\\l 2 ] 9 ^)\\{,V - VT)VT\\H^gn) 

j: l|WrHlf.L(r»IIVrHlit.(r)) 

- ||/i^'^^Vrt’||i,2(^^(r)), 

where the last estimate hinges on K-shape regularity of Th (cf. Lemma 2.7, (i)). □ 

Let V G H^{r). For each T G Th, let Vt be the constant from Lemma 5.1. For each T E T 
we dehne the near-held and the far-held part of the double-layer potential Uj^ = by 

“jCt - ^T)hr) and := - Vt) {I - Vt)) ■ (5.11) 

Note that {v — VT)r]T £ H^{T) C and (n —nT)(l — Vt) ^ H^{dD.) so that ^ 

H^{U\dfl) are well-dehned. Since = —1 in fl and .^1 = 0 in we have, for every 
T eTh, the identities 

Usi + VT = fo ^ fo (5-12) 

5.2. Inverse estimates for the near-held part • The following proposition provides 
an estimate for the near-held part of the double-layer potential. 

Proposition 5.2 (Near-held bound for .h). Let Wh be a a-admissible weight function. There 
exists a constant Cnear > 0 depending only on dLl, F, the n-shape regularity ofTh, and a such 
that the near-field part satisfies G 77^(F), u'^Tln ^ H^{^), eLnd u^t\u\q E 

H^{U \ Vt) with 

Ten 

< Cnear l^ft || ^2(r) • (5.13) 

Proof. Recall stability (2.10) of : H^{dVl) H^{dVl). Taking into account (5.2) 

and the Poincare-type estimate (5.5), we observe 

(2.10) (5.5) 

||Vr7cr*'*^“r lUhr) ~ \\{v - VT)vT\\HHdn) ^ IIVrw||L 2 (^^(r)). 
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Summation over all T G 7h shows 


lli2(T) ~ ll'^/i^r'y||^2(r)- (5-14) 

Ten 

Next, we use the continuity of ^ ; H^^^{dQ) —?■ H^(U \ dQ) from (2.7) and get 

llL2(f/ynO) + llL2(f/ynO«^t) ~ II ~ '^'r)7T||^l/2(gQ) 

Summation over all T G 7h gives 

||wh/h^/^||ioo(-r)(^11 VM““||i2(^j,nQ) + II^'*^^!t lli2(^j,no'=’'*)) ~ Il'^/i^r'y||i2(r)- (5.15) 

Ten 

Combining (5.14)-(5.15), we conclude the proof. □ 


5.3. Estimates for the far-held part As for the simple-layer potential, we have 

a Caccioppoli inequality for the double-layer potential, which underlies the analysis of the 
far-held contribution. 

For the next result, recall from (5.1). 

Lemma 5.3 (Caccioppoli inequality for constant Ccacc of Lemma f.3, the func¬ 
tions of (5.11) satisfy ^ C°°(h), u^^bext g C°°(f2®^*), and u^^lrb^ G C°^{Uf) 

together with 

\\D^u^nn\\LHBsum/six)) < Ccacc II VmJt lbhs,MT)/ 4 W) all X E T E Th. (5.16) 

Proof. The proof is very similar to that of Lemma 4.3. One observes that the far-held 
solves the transmission problem 

= 0 inOuO®^*, 

Wkn] = ~ 'i^T)(l — hr) in i7^'^^(dO), 

[^lU^^^rp] = 0 in i7“^/^(dO). 

We note that (1 — hT)|rn; 7 ' = 0 by construction of 7]t in (5.2). Hence, the same reasoning 

as in the proof of Lemma 4.3 can be applied to reach the conclusion (5.16). □ 

Lemma 5.4 (Local far-held bound for .h). For allT E Th 

l|fc'-'Wnj“«|'TlU>(T) < l|fc"'W4y|U,(T) < C(„ IIv4“,.|U2,„„, (5.17) 

The constant Ciar > 0 depends only on dQ and the n-shape regularity constant ofTh- 

Proof. The lemma is shown in exactly the same way as the corresponding bound for the 
simple-layer potential 53 in Lemma 4.4, appealing to the Caccioppoli inequality (5.16) instead 
of (4.13). □ 
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Proposition 5.5 (Far-field bound for ^). Let Wh be a a-admissible weight function. There 
is a constant Cfar > 0 depending only on dQ, the K-shape regularity constant of Th, and a 
such that 


11 11 L°° (T) 11 ^11 i2 pp) 

Ten 

Ten 

< C'far (lkhVru||i2(r) + lkh//i^^^lli-(r)ll'«^ll|i/2(r)) • 

Proof. Lemma 5.4 implies 

TGTft TGTft Ten 

~ ll'^/i/^^'^^lli°°(T)ll^'*^tiTllL2([/^noe’'t)- (5.18) 

TeTh TgT,, 

With the identities (5.12) and a triangle inequality, we therefore obtain 

||u’hVM^y||^2(y) 

Ten 

( 5 - 18 ) / ~ ~ \ 

^ lk/xA^^li-(T)(l|Vi?(n-nr)||i2(t;^no) + l|Vi?v||i2(t;^no-‘)) (5.19) 

Ten 

Ten 

The near-held contribution is bounded by Proposition 5.2. Furthermore, noting V^vt = 
V(— ut) = 0 in fl, we get 

Ten 

( 5 . 13 ) / ~ ~ \ 

< \\Wh/h^^YL^{T) Y. ( II V.^n||^2([/^no) + ll^•^^llL2(^7^n!^^™p j + \\'^h^rv\\‘j^2(^p') 

Ten 

(2.7) 

< |kft/h^/2||i,x=(r)|lw|||i/2(r) + lkhVrw||i2(r), 

where we have used continuity (2.7) of .ft, the overlap property (4.4), and ||n||/)^i/ 2 (r) = 
Iklli7i/2(an)- n 


6 . Proof of Theorem 3.1 


Finally, we are in a position to prove the inverse estimates (3.1), (3.2) of Theorem 3.1. 


Proof of the inverse estimate (3.1). Let fj G L^(F), extend by zero to the entire boundary 
dQ, and recall the notation from Section 4.1. First, we treat the simple-layer potential QJ. 

22 


With the bounds of Propositions 4.2 and 4.5 we get 


T&Th 

< Y lk/*Vr7r<Tlli3(T) + E W^^r^rYo'CTllhiT) (6-1) 

T&Th TeTh 

< ||W^^^^llL-(r)IIV’ll|-i/2(r) + W^h^h^r)- 

The estimate for the adjoint double-layer potential follows by similar arguments. We split 
the left-hand side into near-held and far-held contributions to obtain 

\\wh^'^\\‘L‘2{r) 'S ll'^/i||L°°(r)ll-^ (t/’Xc/Tnr)||i2(2^)-f- 'Y^ \\'^f>-\\‘L°°{T)\\^ i'^Xr\UT)\\‘L^{T)- (6-2) 
T£Th T&Th 

The continuity ^ : L‘^{dVL) —)■ L‘^{dQ) stated in (2.11) yields for the near-held contribution 
\\'^h\\‘L°°{T)\\^'{'4’XuTnr)\\‘L^(T) ^ "Y^ ll'^/i|lL°°(T)ll-^^(V'Xc/Tnr)|li2(aQ) 

T&Th TeTh 

~ ll'“^/illL°°(T)llV'|lL2(i7Tnr) ~ ll'^ftV'lli2(r)- 
Ten 

For the far-held contribution, we write = 2J('0xr\c/T) note that —1/2 -|- 

and clearly {'ipXr\UT)\T = 0. Therefore, on T we have ^{'iIjxv\Ut) — Furthermore, 

by the smoothness of near T (see Lemma 4.3), we have on T (cf. 

Remark 2.2) and get 

\\^{'fpXT\UT)\\L^(T) = WlT'^^rWh'^iT) = \\duU^nT\\L^(T) ^ 11 I(T’) • 

The far-held contribution in (6.2) can therefore be bounded by Proposition 4.5 

'Yj (V'Xr\nT)llL2('r) ^ 'Y^ \\'^hSu^^p\\\2{p'^ 

Ten Ten 

^ WwhijWhir) + 

Altogether, this gives 

Wwh^'tpWniv) < Wwh'tpWniv) + ||wh/h^/^||L-(r)||^/|l5-i/2(r9 

Proof of inverse estimate (3.2). First, we treat the double-layer potential Let v G iL^(F), 
extend v by zero to n G H^{dQ), and recall the notation from Section 5.1. We recall the 
stability of .^ = ^ : I{^(df2) —)■ II^(df2), from which we conclude G iL^(F). 

Therefore, 

||'M^ftVr^w|lL2(r) = ||rChVr(| -|-7“’'.ft)n||L2(r) < -||'?i’/iVrt'||L2(r) -|- ||w/iVr7o“*Mj?||L2(r) (6-3) 
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with = ^v. There holds Uj^ + vt = in f2, cf. (5.12). For the second term on 

the right-hand side in (6.3), we obtain 

\\whVnTusi\\l2^r) < II^^A^^^IIi-(T)II^^^^Vr7r(“J? + '«^r)||i 2 (r) 

T€Th 

(5.12) 

~ ll'^/i/^^'^^llL°°(T)ll^^^^^r7cr*'*^^!r IIl 2 (t) (6.4) 

Ten 

Ten 

The hrst sum can be bounded by Proposition 5.2, whereas the second sum can be bounded 
by Proposition 5.5. Altogether, this yields 

||r(;/iVr.^n||i 2 (r) < \\wh/h^^‘^\\L^{v)\\'v\\fji/2^Y) + ll'^/iVrw||i, 2 (r) 

and concludes the hrst part of the proof. 

The result for the hypersingular integral operator 211 is shown with similar arguments. Let 
again v G H^(r) and vt as in Lemma 5.1. Note that 2Iln'r = 0. Splitting now into near-held 
and far-held yields 

\\whWv\\l2(^r) = Y1 ll^^2n(n-nT)||i2(T) 

Ten 

< \\whW{{v - vt)vt)\\12(^t) + Y ll^^2n((n - VT)il - VT))\\l2^Ty (6-5) 
Ten Ten 

The near-held contribution is bounded by the stability of QU ; H^{dQ) —)■ L‘^{dQ) stated in 
(2.12) and the Poincare-type estimate (5.5) 

(2.12) (5.5) 

|| 2 IJ((n-nT)hT)|li 2 (r) < \\{v-vt)vt\\Ii(^^(^t)) ^ IIVr^^llL 2 (^,(T))- 

The sum over all elements gives 

||u'/iQIl((n — n'r)hT)Ili2('r) < 'Y^ ll'^^llL°°(T)ll^r'y|li2(^^(2’)) ^ ll'^feVr'y||i2(r)- 
Ten Ten 

It remains to bound the second term on the right-hand side in (6.5). In view of the support 
properties of rjT, the potential = ^{{v — nr)(l — tjt)) is smooth near T (cf. Lemma 5.3) 
so that ~ Furthermore, since QU = — 7 “*.^ we see 

||2U((n — nT)(l ~ 7 T))||L 2 (r) ~ ll 7 r*^'^tirllL 2 (T) = ll^^^'^^TllI, 2 (T) ^ ll^'*^^,Tlli 2 ('r)- 
We use Proposition 5.5 to conclude 

||n;,2U((n-nT)(l-7T))||i2(r) < ||'?r'/iVM^y||^2(7’) 

Ten Ten 

^ lkhVrn||i2(r) + |kh/h^/licx>(r)||n|||i/2(r)- 

Altogether, we obtain 

||M;feQUn||i2(p) < ||wfeVrw||L2(r) + lk/i/^^^^IU“(r)||'y||jpi/2(r). □ 


24 


Acknowledgement. The authors are supported by the Austrian Science Fund (FWF) 
under grant P27005 (MF, DP) as well as through the FWF doctoral school W1245 (MF, 
JMM, DP). The are also supported by CONICYT through FONDECYT projects 3150012 
(TF) and 3140614 (MK). 


Appendix A 


Lemma A.l. Let Th be a regular, n-shape regular triangulation o/F. Suppose that d>2 
and that is a a-admissible polynomial degree distribution with respeet to Th- Then, there 
exists a eonstant Cinv > 0, which depends solely on dQ, the n-shape regularity ofTh, and a, 
such that 

\\h^/\qh + 1)-' ^/^||L 2 (r) < a,v||T,||^-v 2 (r) for all e V^iT)- (A.l) 


Proof. For each T G Th, let XT,qu{T) be a smooth function on with the following properties 
for some hxed 5 > 0 (see, e.g., the proofs of [GeoOS, Lem. 3.7, Prop. 3.8] or the arguments 


below): 

suppxT,q;,(r) Y {x e Tref: dist(a;, 9Tref) > 6/{qh{T) + 1)^}, (A. 2) 

0 < Xr,g;,(T) < 1 in Tref, \\'^XT,qh{T)\\L°°{T,,{) ^ {qh{T) (A.3) 

XT,qUT) = 1 in {x e Tref: dist(a;, ^Tref) > 3S/{qh{T) + 1)^}, (A.4) 

h\\L^T,,() < C\\TTXT,p{T)\\L^T,,t) V polynomials vr of degree qh{T), (A.5) 


lkxr,g^(T)||rfi(Tref) < <^(1 + Qh{T)f\\TT\\L^(T,,i) V polynomials tt of degree g/r(T). (A.6) 


XT,qh{T) is obtained from a mollihcation of the characteristic function of Tref \ 'S' 25 /(g^(r)+i) 2 , 
where := {x G Tref: dist(a;, 9Tref) < e}. The parameter 5 > 0 is dictated by the re¬ 
quirement (A.5). For this, we use the shorthand e{6) = 36/{qh{T) 1)^ and observe 

that we assume XT,qh{T) = 1 on Tref \ S^(^s) so that we are done once we have established 
IKIIl 2 (s^( 5 )) ^ IKIUhnref\Se( 5 ))- [LMWZIO, Lemma 2.1] and the polynomial inverse estimate 
Ikllrfi(Tref) ^ iQh{T) + l)2]|7r]|L2(T^^j), yield 

Mlnseis)) ~ < £{S){qh{T) + 


= emqh{T) + ir 



2 

nhTef\5,(«)) 





Taking 6 sufficiently small produces ]|7r]|f,2(5^^^^) < l|7r]|L2(Tref\Se(5)) ns desired. 

Dehne XT,qh{T) with suppxr, 5 ^(T) ^ T by Xr,q^(T) o 7 t = XT,qh{T)- Given ^h e T^iT), dehne 
H^(T) 3 V := 'YfjT&Th elementwise fashion by requiring suppu^ Y T and 


Vt\t 


h{T) 

{^ + qh{T)y 


h\T)XT,qH{T), 


(A.7) 


Note that Vt G H^fT) by the support properties of XT,qh(T)- An interpolation inequality and 
the estimate (A.6) on the reference element give 


ll^rll|i/ 2 (r) 


< Ibrll 


|i/2(9q) ^ llwTl|L2(an)llwTl|Hi(aO) - \\vT\\L2{T)\\vT\\m{T) 


(A.6WA.5) (1 + g^ 

~ h(T) 
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This implies 


2 


T&Th 

From [SSll, Lemma 4.1.49] we get from (A.8) 

2 


+ 1 


hV2 


L2(r) 


iHi/2(r) 


Finally, we estimate 


T€Th 


< 


E ii-tIi 


, < 

rri/2(r) ~ 




1 + g/i 


hV2 


L2(r) 


(A.8) 


(A.9) 


^1/2 2 


1 + g/i 


= E 

L2(r) 


/i(T)1A 


1 + (lh{T) 




(A.5) 

L^{T) rp&Th 


h(T)iA 2 

-XT,qh{T)^h 


1 + (lh{T)' 


L^{T) 


(A.2) 

< 


h(T)V2 




^\/XT9^^^lli2(T) - (^^T, ^/i)l2(t) - (n, T/i)i2(r) 

^ ’ T&Th 


(A.9) 


— ll^fcll/r-i/2(r)||'i^|lHi/2(r) ^ ll^/i||H-i/2(r) 


1 + g/i 


hV2 


L2(r) 


(A.7) 

< 11^ 


/j|li7-i/2(r) 




1 + g/i 




□ 


L2(r) 
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